Multispecies Ion Acceleration in 3D Magnetic Reconnection with Hybrid-Kinetic Simulations.

Magnetic reconnection drives multispecies particle acceleration broadly in space and astrophysics. We perform the first 3D hybrid simulations (fluid electrons, kinetic ions) that contain sufficient scale separation to produce nonthermal heavy-ion acceleration, with fragmented flux ropes critical for accelerating all species. We demonstrate the acceleration of all ion species (up to Fe) into power-law spectra with similar indices, by a common Fermi acceleration mechanism. The upstream ion velocities influence the first Fermi reflection for injection. The subsequent onsets of Fermi acceleration are delayed for ions with lower charge-mass ratios (Q/M), until growing flux ropes magnetize them. This leads to a species-dependent maximum energy/nucleon ∝(Q/M)^{α}. These findings are consistent with in situ observations in reconnection regions, suggesting Fermi acceleration as the dominant multispecies ion acceleration mechanism.

Efficient Nonthermal Ion and Electron Acceleration Enabled by the Flux-Rope Kink Instability in 3D Nonrelativistic Magnetic Reconnection.

The relaxation of field-line tension during magnetic reconnection gives rise to a universal Fermi acceleration process involving the curvature drift of particles. However, the efficiency of this mechanism is limited by the trapping of energetic particles within flux ropes. Using 3D fully kinetic simulations, we demonstrate that the flux-rope kink instability leads to strong field-line chaos in weak-guide-field regimes where the Fermi mechanism is most efficient, thus allowing particles to transport out of flux ropes and undergo further acceleration. As a consequence, both ions and electrons develop clear power-law energy spectra that contain a significant fraction of the released energy. The low-energy bounds are determined by the injection physics, while the high-energy cutoffs are limited only by the system size. These results have strong relevance to observations of nonthermal particle acceleration in space and astrophysics.

Scaling of magnetic reconnection in relativistic collisionless pair plasmas.

Using fully kinetic simulations, we study the scaling of the inflow speed of collisionless magnetic reconnection in electron-positron plasmas from the nonrelativistic to ultrarelativistic limit. In the antiparallel configuration, the inflow speed increases with the upstream magnetization parameter σ and approaches the speed of light when σ>O(100), leading to an enhanced reconnection rate. In all regimes, the divergence of the pressure tensor is the dominant term responsible for breaking the frozen-in condition at the x line. The observed scaling agrees well with a simple model that accounts for the Lorentz contraction of the plasma passing through the diffusion region. The results demonstrate that the aspect ratio of the diffusion region, modified by the compression factor of proper density, remains ∼0.1 in both the nonrelativistic and relativistic limits.

Formation of hard power laws in the energetic particle spectra resulting from relativistic magnetic reconnection.

Using fully kinetic simulations, we demonstrate that magnetic reconnection in relativistic plasmas is highly efficient at accelerating particles through a first-order Fermi process resulting from the curvature drift of particles in the direction of the electric field induced by the relativistic flows. This mechanism gives rise to the formation of hard power-law spectra in parameter regimes where the energy density in the reconnecting field exceeds the rest mass energy density σ ≡ B(2)/(4πnm(e)c(2))>1 and when the system size is sufficiently large. In the limit σ ≫ 1, the spectral index approaches p = 1 and most of the available energy is converted into nonthermal particles. A simple analytic model is proposed which explains these key features and predicts a general condition under which hard power-law spectra will be generated from magnetic reconnection.

Nonlinear evolution of the lower-hybrid drift instability in a current sheet.

The lower-hybrid drift instability is simulated in an ion-scale current sheet using a fully kinetic approach with values of the ion to electron mass ratio up to m(i)/m(e)=1836. Although the instability is localized on the edge of the layer, the nonlinear development increases the electron flow velocity in the central region resulting in a strong bifurcation of the current density and significant anisotropic heating of the electrons. This dramatically enhances the collisionless tearing mode and may lead to the rapid onset of magnetic reconnection for current sheets near the critical scale.